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a— decay through a deformed potential barrier produces significant mixing of angular momenta 
when mapped from the nuclear interior to the outside. Using experimental branching ratios and 
either semi-classical or coupled-channels transmission matrices, we have found that there is a set 
of internal amplitudes which are essentially constant for all even-even actinide nuclei. These same 
amplitudes also give good results for the known anisotropic a— particle emission of the favored decays 
of odd nuclei in the same mass region. 
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, The phenomenon of "tunneling in the presence of an environment" is of considerable interest in many branches 
iy~j ■ of physics and chemistry |Q, for example in the tunneling of Cooper pairs through a Josephson junction [0. In the 
domain of nuclear physics, the study of the sub-barrier fusion of heavy nuclei has made significant contributions to 
^sO ■ this problem over recent years, in that the experimental distribution of fusion barriers has been shown to be intimately 
related to the environment consisting of the target and projectile excited states S. In this problem the incident wave 
is known and all the transmitted flux ends up in the single fusion channel. The phenomenon of a— decay is potentially 
more difficult, since the incident wave (a— particle wave function in the nuclear interior) is unknown but potentially 
much more rewarding since (a) for even-even nuclei the transmitted flux may end up in different daughter states 
for which the individual fluxes, i.e. branching ratios, can be measured and (b) since for odd nuclei, the presence of 
different angular momenta in each daughter state may lead to a measurable anisotropy in the a— emission. 

The a— decay of deformed nuclei may be divided into two distinct parts: the formation of an a— particle in the 
. £h ! nuclear interior, followed by its penetration through the a— daughter deformed Coulomb barrier. There are various 
approaches to the formation problem. One of these assumes a preformed a— particle (or at least a spatially correlated 
4-nucleon cluster with the appropriate quantum numbers), which moves in the deformed field of the daughter nucleus 
H,[|. Froman |6| assumes a constant a— particle probability on the deformed nuclear surface. Mang et al. [Q and, 
more recently, Delion et al. ^ have considered the deformed single-particle (Nilsson) states in the vicinity of the 
Fermi surface and taken the overlap of the correlated neutron and proton BCS wave functions with an a— particle at 
the nuclear surface. 

Whatever the formation mechanism, the decay proceeds by penetration through the Coulomb barrier. If this barrier 
is assumed to be spherical, there will be no mixing of orbital angular momentum states during the tunneling and, 
for an even-even nucleus, the observed branching ratios to different rotational states of the daughter (/ = + , 2 + , 4 + , 
....) will be determined by the L— admixture in the nuclear interior, modified by the transmission factors for the 
different centrifugal barriers evaluated at slightly different a— particle energies due to the excitation energies e} of the 
daughter. If one takes account of the deformation of the barrier, there will be additional mixing during the tunneling 



Such effects have been considered in the fusion of heavy nuclei |10| | and more recently confirmed experimentally 
[U|,|l2|. In this Letter, we take known branching ratios and calculate the mixing in the barrier to obtain the internal 
amplitudes near the nuclear surface. No model of internal dynamics (a— preformation or distribution) is required 
and indeed one does not know what preformation factors have to be fitted unless the barrier penetration problem is 
first solved. We then systematically survey the internal amplitudes for the even-even actinide a— emitters for various 
choices of relative phases. A surprising result that emerges from our analysis is that the relative internal amplitudes 
may possibly be constant over a wide range of actinide nuclei and this in itself gives a very strong indication of the 



1 



type of nuclear model for the preformation factor which might be the most appropriate. 

The fine structure of the a— particle energy spectrum determines the branching ratios to states of known spins I 
in the daughter. Since, for an even-even nucleus, the orbital angular momentum of the a— particle must be L = I, 
we shall set up our problem in terms of these outgoing L— waves. In the intrinsic frame of the axially-symmetric 
deformed daughter, one may write an asymptotic outgoing wave in the form 

r utslde =J2 C LO° L utslde (k L ,v) = J2°l {GL(kL,r)+iF L (k L ,r))Y LO (r)xLo, (1) 

L L 

where Fl and Gl are the regular and irregular Coulomb wave functions. The Q L utmde represent outgoing a— particles 
with orbital angular momentum L coupled to total angular momentum zero with a state xlo of the daughter of 
spin I = L. The wave numbers fcj, differ due to the different e*j. With such an outgoing wave, one can perform 
a coupled-channels calculation to obtain the wave function on the other side of the barrier. This must have both 
incoming and outgoing components 

^inside = J2 (A L Of side (k L ,r) + B L 2 L nside (k L ,r)) . (2) 

L 

The corresponding currents are shown schematically in Fig. 1. 
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FIG. 1. Schematic diagram of the potential V and the incoming and outgoing currents. The classical turning points n and 
r2, and the energy of the a— particle E a are used in the semi-classical calculation of the transmission matrix, see Eq. (JE]). 

The above problem may be solved numerically to obtain the outside coefficients Cl in terms of those inside, Al- 
One may thus define a transmission matrix M 



C L =^2M LL , A l , 



(3) 



Since we wish to undertake a systematic study of a— emitters in the actinide region, we shall first obtain the trans- 
mission matrix using the semi-classical method outlined below. We thus express the outer amplitudes ci (lower case 
indicates WKB solutions) in terms of the amplitudes in the nuclear interior by 



cl = y^JC LL <a L ' = y^(Y L0 1 exp(- I L (cos 0))\Y L > o )a L r 

L> L' 

where I takes the usual WKB form 

Il{cos6) -- 



n (cos o) 



2/t 

Ti 



2 (V L (r,cos6)-E a ) 



dr . 



(4) 



(5) 
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(see Fig. 1). The potential Vl(t, cos 6) comprises the Coulomb field due to the deformed charge distribution of 
the daughter, a deformed Woods-Saxon potential and a centrifugal term. We consider both fa and fa terms in the 
deformation. The angular integrals were evaluated using the technique of Kermode and Rowley |l3| ]. The matrix 
K-ll' is the analogue of Mll> m the coupled-channels formalism. 

The above approach is essentially an exact calculation of the semi-classical transmission coefficients introduced by 
Froman || and used, for example, in Ref. ||. Indeed the earlier expressions of Froman may be obtained from Eq. (Q) 
by ignoring the hexadecapole deformation, taking the nuclear potential to be a 9— dependent sharp cut-off and by 
making a first-order expansion in fa . Approximations similar to these have also been employed in a coupled-channels 
formalism 

The magnitudes of the coefficients are proportional to the square roots of the branching ratios for the angular 
momenta L. They may be taken to be real and can, in principle, be cither positive or negative. Their reality follows 
from the requirement that the imaginary part of the wave function at the nuclear surface should be small [Q. In 
this Letter, we restrict our calculations to {L} = {0,2,4} since few nuclei have measured branches to the 6 + state. 
(We have, however, considered the inclusion of 6 + and 8 + states and have found that our conclusions are essentially 

unaffected.) We consider 4 possibilities for the relative signs of the {cl}, i.e. {+ + +}, {+H — }, {H h} and {H }. 

For all the decays we consider, the Sommerfeld parameters are large and the Coulomb phases then ensure that for 
the case {+ + +} the spherical harmonics in Eq. ([!]) are in phase along the symmetry axis. Thus if one considered 
the amplitudes in Eq. (Q) to add coherently, the outgoing flux would be axial. In the even-even system, of course, 
the currents are not added coherently and the outgoing flux is always isotropic for each L, after integrating over all 
orientations of the daughter. However, we shall sec below that the above consideration is important for odd-even 
systems which may be polarized to yield anisotropic a— decay. For each of the above phase combinations, we have 
determined the amplitudes from Eq. (|4|) using deformation parameters from Moller et al. p6| . Figure 2 shows 
the 4 sets of (normalised to unity) for actinide nuclei with daughter mass (atomic number) 220 < A c i < 248 
(88 <Z d < 96). 
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FIG. 2. The amplitudes a_L for even-even actinide nuclei with the four choices of external phase indicated. The circles, 
squares and triangles represent L = 0, L — 2 and L — 4, respectively. 



To test the accuracy of the WKB approximation, we have performed exact coupled-channels calculations for three 
of the above nuclei. The results were in good agreement with those shown in this figure. For example, for 238 U 

with the phase choice {cl} = {H }, the coupled-channels give {a^} = {0.84,-0.54,0.08} compared with the 

semi-classical values {ol} = {0.83, —0.55, 0.10}. The importance of the L— mixing under the barrier is demonstrated 
by the corresponding results {a^} = {0.70, —0.68, —0.19} for the spherical case. 

Since the deformation of the daughter varies smoothly with Ad, one might expect that just one set of the solutions 
shown in Fig. 2 corresponds to the physical amplitudes. One particularly interesting solution is that obtained with 

the combination {H } (Fig. 2d) since the coefficients are practically nucleus-independent even though fa varies 

from around 0.10 to 0.24 and the a— particle energies vary from around 4 to 7 MeV over this mass region. It is also 
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the only one to yield amplitudes with |ao| > \ai\ > for all nuclei (see the discussion below). 

We note, however, that for the constant solution <22 is always negative. At first sight, this implies that if the same 
solution prevailed for odd-even nuclei, then a— emission takes place mainly in the equatorial plane - a result which 
would not be consistent with measured anisotropics. Consider, however, the favored a— decay of an aligned nucleus 
with K = M = J p = Jd, where K is the spin projection of the nucleus along the symmetry axis, M the projection of 
the parent spin J p along the space-fixed z— axis and Jd the daughter spin. For this case the angular distribution of 
the emitted a— particles may be written (see, for example, |^Jl^|) 

W{9) oc I J2( Jd LM - m m \J P M)(J P L K 0\J d K) c L F im (0,O)| , (6) 

m L 

where the Clebsch-Gordon coefficients arise out of the transformation from the body-fixed to the lab frame. We define 
the anisotropy by W(0)/W(^tt). Equation (||) describes the favored decay K p = = K angular distribution only. 
A contribution from the decay K p = K — > K d = —K is suppressed (see, for example, page 50 of |J or pages 272-3 of 

A simple model for the odd-even nucleus would be to assume that it consists of an even-even core plus a spectator 
nucleon. Indeed Froman || employed the above formula using cl extracted from the neighbouring even-even branching 
ratios, which implicitly uses the same transformation (^) as for the even-even case. Delion et al. [0 have also used 
Eq. (|J) but derive their amplitudes cll from a BCS calculation. The expression (fj) may not, however, be appropriate 
for odd-even systems, where the a— particle energy is determined by the daughter spin Jd rather than by the orbital 
angular momentum L. We have applied Eq. (Q) to the four odd-even actinide nuclei for which anisotropies have been 
measured (see Table |), using the amplitudes from Fig. 2. We find that Eq. (||) cannot predict anisotropies both 
greater than and less than one for these nuclei, for any of the four sets of Of,. This is because the sign of cijc^ is 
either positive (anisotropy > 1, i.e. Figs. 2a, 2b) or negative (anisotropy < 1, i.e. Figs. 2c, 2d). 

If the energies of the excited states of the daughter were high, then barrier penetration would filter out the compo- 
nents of the wave function corresponding to configurations other than its ground state. In that case we obtain 

c L ~^K-ll-{J v V K 0\J d K) 2 a L ,. (7) 
v 

We note that Bohr and Mottelson Q suggest a similar J^-decoupled equation for the leading order transition rates. 
Equation (Q) has also been applied to the four nuclei in Table Q, using the internal amplitudes of Fig. 2d. 

TABLE I. Anisotropies {W{0) /W '(5^)), for the favored decay of four odd -even nuclei. Deformation parameters correspond 
to the daughter nucleus. 

Parent J p @2 Pa Theory {H } Experiment ref. 



221 j 

227j 

241 Am 

253 Eg «* 



Fr 
7 Pa 



0.039 
0.147 
0.215 
0.235 



0.028 
0.110 
0.102 
0.040 



0.77 
2.66 
4.26 
3.70 



0.37(2) 
3.55(28) 
>2.7 
>3.8 



1121 
1121 

[HI 



Moller et al. give a non-zero value for j3s 
* includes L=6 



The values of the predicted anisotropies agree well with the experimental data. In particular, we emphasize the 
prediction of anisotropies both less than and greater than 1, i.e. it is possible to have more a— particles emitted along 
the symmetry axis than equatorially, despite the fact that 02 < 0. The reason is that the external amplitude C2 may 
become positive, since the Clebsch-Gordan coefficients of Eq. (fj]) attenuate the effect of the internal amplitude 02. 
We note that this effect is not possible if we use the amplitudes from Figs. 2a or 2b. The solutions of Poggenburg 
et al. and of Delion et al. give anisotropies always greater than one. However, the excited states of the daughter lie 
at relatively low excitation energies and, in the present model, their coupling to the ground state is not sufficiently 
attenuated by barrier penetration for the reversal of sign between ai and C2 to take place. The success of Eq. (|t]), 
could, however, suggest some other dynamical element through which the a— particle orbital angular momenta are 
mixed in exactly the same way as for even-even nuclei, but the coupling to different daughter states is absent. 
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We have shown that it is possible to describe all known branching ratios of even-even actinide nuclei with an 
a— particle wave function near the nuclear surface which is practically nucleus-independent. This model has a certain 
aesthetic appeal in itself and, with the assumption that the same L-mixing matrix JCll 1 is present in the favored 
decay of an odd system, is capable of reproducing the known anisotropies in this mass region. This of course begs 
the question as to what physical model could generate such constant amplitudes. The best candidate would appear 
to be the notion that the a— particle amplitudes should be projected from the pair-correlated neutron and proton 
Nilsson-model states. In this mass region, the level density is high and the pair forces lead to a rather diffuse Fermi 
surface. One might then expect that the correlated ground-state wave function should vary rather slowly with the 
Fermi energy, or in other words with the nucleon numbers of the system. Since pairing mainly takes place through 
the two-body angular momentum J = 0, this model would also be expected to give amplitudes with the property 
Oo| > a 2 > |a-4 1 , as found in Fig. 2d. 
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